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1  Introduction 


The  use  of  feedback  controllers  that  satisfy  the  positive  real  property  are  known  to 
have  robust  stability  properties  and  have  been  shown  to  be  useful  especially  in  regard 
to  the  control  of  flexible  spacecraft  [1,  2,  3].  In  this  case  the  positive  real  controls 
can  prevent  destabilization  of  the  closed  loop  system  due  to  control  or  observation 
spillover,  and  due  to  erroneous  modal  information.  In  a  previous  paper  [3]  we  applied 
a  positive  real  synthesis  technique  to  the  design  of  large  space  structure  controls.  In 
this  case  the  controller  eigenstructure  was  determined  by  ad  hoc  manipulations  of  an 
arbitrary  symmetric  positive  definite  matrix.  To  obtain  controller  eigenvalues  with 
“sufficient”  damping  to  give  a  reasonable  closed-loop  response  was  difficult  in  some 
cases  when  using  a  diagonal  form  for  this  matrix.  A  more  general  arbitrary  matrix  is 
desired. 

In  this  report  we  discuss  a  technique  for  the  positive  real  controller  design  that 
allows  general  eigenvalue  placement,  hence  resulting  in  improved  design  flexibility. 
This  technique  is  based  on  recasting  the  eigenvalue  placement  problem  as  a  con¬ 
strained  optimization  problem.  Two  formulations  for  the  solution  of  this  problem  are 
developed,  and  numerical  results  are  presented  for  both  algorithms. 


2  Positive  Real  Systems 

A  time  invariant  single-input  single-output  system  is  said  to  be  positive  real  if  the 
scalar  system  transfer  function  G{s )  satisfies  the  inequality  3?[G(s)]  >  0  for  all  -R(s)  > 
0,  and  the  system  is  strictly  positive  real  if  3?[G(,s)]  >  0  for  all  3?(s)  >  0.  For  a 
multivariable  system,  the  Hermitian  component  of  the  frequency  response  matrix 
must  be  positive  semi-definite  and  positive  definite,  respectively.  For  a  scalar  system, 
positive  realness  implies  the  Nyquist  plot  lies  strictly  in  the  first  and  fourth  quadrants. 
It  is  known  that  a  positive  real  system  is  “dissipative”,  that  is  (generalized)  energy  in 
such  a  system  is  bounded  above  by  the  initial  energy  plus  the  control  energy  [4].  It 
is  this  dissipative  nature  of  the  positive  real  controller  which  makes  it  an  attractive 
choice  for  the  control  of  a  conservative  lightly  damped  structure. 


The  important  stability  theorem  used  to  guarantee  stability  was  first  proposed  by 
Popov  [5].  The  stability  results  can  be  stated  as: 

Popov  Stability  Theorem 

A  sijsh  in  with  forward  transfer  matrix  <  » ( s )  and  fudlmck  matrix  //(s)  fs  (junr- 
autre.d  stable  if  one  of  the  matrices  is  strictly  positive  real  and  the  other  is  positive 
real. 

For  a  flexible  structure  the  transfer  matrix  between  ideal  colocated  actuators  and 
velocity  sensors  is  known  to  be  positive  real,  hence  a  strictly  positive  real  feedback  will 
lead  to  a  stable  closed  loop  response.  While  realistic  actuators  and  velocity  sensors  will 
not  in  general  have  a  positive  real  transfer  function  matrix,  Slater  et  al.  [6]  have  shown 
how  actuator  dynamics  can  be  viewed  as  a  multiplicative  perturbation  to  a  baseline 
positive  real  model.  These  uncertainties  can  then  be  treated  using  the  robustness 
results  of  Doyle  et  al.  [7,  S,  9],  to  derive  sufficient  conditions  for  guaranteed  closed- 
loop  stability.  This  stability  will  then  hold  even  though  the  plant  transfer  function 
matrix  will  not  be  positive  real  over  the  entire  frequency  spectrum. 

2.1  Control  Design 

The  design  of  the  positive  real  feedback  controller  is  described  in  McLaren  [3]  and 
Slater  [10].  Summarized  briefly  here,  the  controlled  system  is  assumed  to  be  in  the 
form 


x  =  Fx  +  Gu 


(1) 


y  =  Hx 


(2) 


A  positive  real  controller  is  designed  using  an  observer  model  based  on  a  reduced- 
order  model  of  the  controlled  system.  The  reduced-order  system  is  assumed  to  be  in 
the  form 


i  r* 


(3) 


2 


(4) 


V  m  X  m 


where  Fm,  Cm,  and  Hm  are  the  matrices  of  the  reduced-order  model,  xm  is  the 
n-dimensional  reduced  state  associated  with  the  reduced-order  model,  y  is  the  p- 
dimensional  output  vector,  and  u  is  the  p-dimensional  input  vector.  The  dynamics 
of  the  estimator  are  given  by 

x  =  Fmx  +  Gmu  +  I\F(y-Hmx)  (5) 


while  the  feedback  control  function  is  assumed  to  be  of  the  form 

u  =  —Kx  (6) 

The  feedback  gain  matrix  K  can  be  chosen  by  any  desired  means.  For  our  applica¬ 
tions,  K  is  chosen  to  minimize  an  index  of  performance 

J  (yT Ay  +  uT Bu^j  dt  (7) 

where  the  weighting  matrices  A  and  B  are  selected  so  as  to  place  the  closed-loop 
poles  of  the  reduced  model  (3)  and  (4)  at  acceptable  locations  (i.e.  to  place  the  poles 
of  [Fm-GmI<]). 

The  observer  transfer  matrix  relating  the  feedback  variable  u'  =  —u  to  the  output 
y  can  be  found  from  equations  (5)  and  (6)  as 

He(s)  =  K(sI-Ft)~lKr  (8) 


with 


Ft  =  (Fm  -  GmK  -  KFHm) 


(9) 
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The  matrix  I\p  is  an  (ti  x  p )  matrix  of  estimator  gains.  He(. s)  represents  the  com¬ 
pensation  transfer  matrix  and  must  be  chosen  to  yield  closed-loop  stability  for  the 
actual  system  (1  and  (2).  He(s )  is  guaranteed  to  be  positive  real  if  Kp  is  chosen  to 
satisfy  [10] 


PFt  +  F?P  =  -Q  (10) 

I\p  =  P~x  I\T  (11) 


with  both  Q  and  P  positive  definite  symmetric.  Equation  (10)  may  be  rewritten  as 
PFd  +  FjtP  =  -  (Q  -  I<THm  -  HlK)  =  -Q  (12) 


where 


Fci  =  (Fm  -  GmI\) 


(13) 


is  now  assumed  to  be  known.  Calculation  of  I\p  given  I\  and  Q  follows  directly  from 
equations  (11,  (12)  and  (13). 

The  well  known  properties  of  the  Lyapunov  equation  (10)  specify  that  with  posi¬ 
tive  definite  Q  and  P  we  are  guaranteed  that  Ft  is  a  stability  matrix.  In  our  solution 
technique  using  equations  (6)  and  (7),  the  choice  of  feedback  gain  matrix  K  leads  to 
an  asymptotically  stable  Fd.  The  choice  of  either  a  positive  definite  Q ,  or  a  positive 
definite  Q  docs  not  guarantee  the  positive  real  condition.  Q  >  0  does  not  guarantee 
that  the  solution  P  of  the  Lyapunov  equation  is  positive  definite  since  Q  may  not 
be  positive  definite.  Similarly,  while  Q  >  0  is  sufficient  for  P  >  0,  this  does  not 
automatically  guarantee  Q  >  0  since  the  matrix  (/\'T Hm  +  A')  will  generally  be 

indefinite.  In  a  solution  technique  the  important  necessary  conditions  to  be  enforced 
are  Q  >  0  and  P  >  0.  For  a  given  choice  of  positive  definite  Q  then,  we  need  to  be 
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mindful  that  the  definiteness  of  Q  and  of  P  must  be  ascertained  separately.  Clearly 
though,  for  “large”  Q,  the  definiteness  conditions  on  Q  and  P  are  easily  met. 

A  choice  of  the  matrix  Q  determines  the  estimator  gain  matrix  AV-  Most  of  our 
previous  studies  have  restricted  this  choice  to  a  diagonal  matrix,  or  at  most  one  with 
ad  hoc  variations  of  a  few  matrix  elements.  It  is  clear  however  that  this  severely  limits 
the  eigenplacement  of  the  controller  and  may  lead  to  unsatisfactory  control  properties. 
Figure  1  shows  the  locus  of  eigenvalues  of  the  nominal  closed-loop  estimator  matrix 
Ft  —  [Fm  —  A' p Hm )  for  a  6</l  order  system  when  Q  is  chosen  as  Q  =  ql,  and  illustrates 
the  restrictions  on  eigenvalue  placement  when  using  a  diagonal  Q. 

A  more  general  technique  is  required  in  order  to  determine  a  rational  approach 
to  the  determination  of  the  Q  matrix,  which  can  lead  to  more  general  eigenvalue 
placement  and  hence  a  satisfactory  control  design.  In  the  present  work  we  give  de¬ 
tails  of  a  numerical  procedure  to  determine  Q  on  the  basis  of  eigenvalue  placement. 
Prescribing  a  full  Q  matrix  involves  choosing  the  n(n  +  l)/2  independent  elements, 
since  Q  is  symmetric.  If  these  values  are  to  be  chosen  to  fulfill  an  eigenvalue  require¬ 
ment,  there  are  only  n  equations  specifying  these  elements,  so  the  problem  is  highly 
under-determined. 

There  are  no  obvious  correlations  between  the  closed-loop  controller  eigenvalues  of 
Fe  and  the  choice  of  Q.  In  general,  the  AV  required  to  achieve  a  given  set  of  desired 
eigenvalues  is  not  unique,  while  the  choice  of  Q  to  achieve  a  given  AV  is  also  not 
unique.  Finally,  it  is  easy  to  show  that  many  eigenvalue  placements  are  not  possible 
within  the  positive  real  constraint. 

3  Pole  Placement  for  the  Positive  Real  Estimator 

The  essence  of  this  work  is  to  show  that  the  non-linear  pole  placement  problem  for 
the  positive  real  estimator  (10),  (11)  can  be  converted  to  an  equivalent  optimization 
problem  for  which  there  are  a  variety  of  numerical  algorithms.  Two  algorithms  are 
presented  which  differ  in  their  choice  of  performance  index  and  in  the  method  of 
handling  constraints. 
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3.1  Eigenvalue  Sensitivity  Method 


Consider  the  case  where  a  nominal  choice  of  Q  in  (12)  results  in  a  given  set  of 
eigenvalues  for  Fe,  which  are  written  in  vector  form  as  A.  These  differ  form  thv 
desired  eigenvalues  A  by  an  amount 


6X  =  A  -  A 


(14) 


where  dirn(A)  =  n.  Let  the  independent  elements  of  Q  be  written  in  vector  form  as 


9  = 


< In  <7i2  <7i3 


T 


(15) 


where  dimfqr)  =  n(n  +  l)/2. 

Small  changes  in  q  will  cause  changes  in  A.  To  achieve  the  desired  eigenvalue 
placement  we  need  to  find  the  solution  6 A  =  f{q),  which  may  be  an  underdetermined 
set  of  non  linear  equations.  If  we  assume  the  changes  in  A  and  q  are  small,  then  we 
can  use  linear  theory  and  write 


<5A  =  M6q 


(16) 


where  ,\I  is  a  matrix  of  partial  derivatives.  M  is  of  dimension  ( n  x  n(n  +  l)/2), 
hence  there  are  a  multiplicity  of  solutions  for  6q.  To  stay  within  the  linear  regime, 
one  approach  is  to  find  the  "smallest”  6q  to  minimize  the  changes  in  Q.  i.e.  the  least 
squares  solution 


6q  =  \I*6\ 


(17) 


where 


A/t  =  A/T(A/A/T)_1 


(IS) 
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[f  the  original  6\  is  large  then  the  6q  given  by  (17)  is  large  and  the  linear  theory 
is  no  longer  valid.  Home  we  must  break  the  actual  6A  into  a  number  of  smaller  steps 
for  which  the  linear  theory  is  applicable,  and  proceed  sequentially  to  the  desired  set 
of  eigenvalues.  Two  problem.*  must  be  solved  to  apply  this  technique.  The  first 
is  to  compute  M ,  the  matrix  of  partials,  and  the  second  is  to  impose  the  positive 
definiteness  constraints  on  P  and  Q  to  ensure  positive  realness  of  the  resulting  transfer 
matrix. 

3.1.1  Eigenvalue  Sensitivities 

The  M  matrix  is  determined  by  evaluating  the  sensitivities  of  the  estimator  eigen¬ 
values  to  parameter  changes.  These  sensitivities  are  found  from  the  properties  of  the 
eigenvalues  and  eigenvectors  of  Fe.  Let  be  the  ith  eigenvalue  of  Fe.  Then  the  right 
eigenvector  u,  and  the  left  eigenvector  v,  satisfy  the  relations 


F,u,  =  A,u,  (19) 

v-F,  =  A,<  (20) 


where  (•)*  indicates  conjugate  transpose,  and  U{  and  v,  are  chosen  to  satisfy  the 
orthogonality  condition 


(21) 

Differentiating  (19)  with  respect  to  the  jkth  element  of  Fe,  premultiplying  by  v‘, 
and  using  (20)  and  the  orthogonality  condition  (21),  we  can  obtain  the  sensitivity 
relationship 


d\, 


(22) 


Using  the  property  that 
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'  OFe 
dFP , 


Sijfikm 


(23) 


we  find  that 


d\,  _  . 

OF  l,jU'k 

Ul  f,  k 


(21) 


where  v"  and  utk  are  the  jth  element  of  vt  and  the  ktfl  element  of  u ,  respectively. 

Using  these  results  we  can  proceed  to  find  the  relationship  of  the  eigenvalue  A,  to  a 
change  in  the  jkth  element  of  Q  through  a  succession  of  chain  rule  operations  involving 
I\  [equation  (11)],  P  [equation  (12)],  and  finally  Q.  Performing  these  manipulations 
we  have  finally 


dXj 

d6(]jk 


v-P 


-i 


dP 

d6qjk 


A  pHm  u , 


(25) 


where  S'jh  is  the  partial  derivative  of  A,  with  respect  to  a  change  in  the  jkth  element 
of  Q,  and 

dP 

d6qjk 

is  a  matrix  obtained  by  differentiating  equation  (12).  This  gives  the  Lyapunov  equa¬ 
tion 


^LFi  +  ft1Z- 

OSqjk  cl  +  c‘  dSq]k 


dQ 

OSqjk 


where  because  of  the  required  symmetry  of  Q 


J^-)  =  {S„S,k)  +  (6rk6„) 

OSqjkJr3 


(26) 


(27) 
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The  matrix  M  in  (16)  is  formed  by  the  appropriate  combination  of  the  sensitivity 
elements  S'jk.  Evaluation  of  the  elements  of  M  as  above  requires  the  solution  of  the 
n(u  +  l)/2  Lyapunov  equations  containe<l  in  (26)  to  obtain  the  sensitivities 

dP 

d6qjk 


Each  equation  is  related  however  since  they  differ  only  on  the  right  hand  side.  Since 
most  techniques  for  the  solution  of  low  order  Lyapunov  equations  require  an  LU 
decomposition,  the  entire  set  of  matrix  sensitivity  solutions  reduce  to  one  LU  decom¬ 
position  and  n(n  +  l)/2  elimination  solutions  to  a  set  of  linear  equations,  thereby 
making  the  numerical  problem  quite  tractable. 

Generally  the  systems  we  deal  with  will  have  complex  eigenvalues,  hence  the  M 
defined  here  is  complex.  To  remain  in  the  real  domain  for  ease  of  computation,  it 
is  convenient  to  redefine  M  and  8X  by  partitioning  (16)  into  real  and  imaginary 
components 


M  = 


(28) 


where  the  subscripts  r  and  i  here  indicate  the  real  and  imaginary  parts  of  the  respec¬ 
tive  quantities.  For  complex  conjugate  pairs  of  poles  only  one  root  is  used,  thus  not 
changing  the  overall  order  of  the  problem,  and  ensuring  that  M  is  not  rank  deficient. 

Also,  if  it  is  not  necessary  to  place  all  n  eigenvalues  of  a  problem,  then  only  those 
eigenvalues  which  are  explicitly  required  to  be  placed  need  be  involved  in  the  analysis. 
Similarly,  if  the  imaginary  part  of  the  estimator  eigenvalues  are  not  required  to  be 
placed,  then  one  can  simply  delete  the  half  of  (28)  that  deals  with  the  imaginary  part 
and  then  solve  the  problem  as  before. 

Finally,  in  the  numerical  implementation  of  this  solution,  care  must  be  taken  to 
ensure  that  the  change  vector  8 A  requested  is  “small",  so  that  one  remains  within 
the  linear  region  (where  changes  in  Q  cause  a  change  in  A  that  is  close  to  the  change 
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predicted  by  the  linear  analysis).  Hence,  in  the  examples  shown  here  we  solve  the 
linear  problem  repeatedly  for  a  number  of  steps  of  the  form  a8 A,  where  the  a  is 
chosen  to  make  the  a6A  small. 

3.1.2  The  Positivity  Constraints 

If,  during  the  eigenvalue  placement  as  already  outlined,  the  Q  matrix  crosses  a  positive 
definite  boundary  and  becomes  and  remains  indefinite,  the  final  He(s)  will  violate  the 
Popov  Stability  Theorem  as  stated  earlier  for  a  positive  real  plant  G(s)  (even  though 
the  matrix  P  may  indeed  still  be  positive  definite!).  Therefore,  at  a  positive  definite 
boundary  of  Q,  we  must  alter  the  change  vector  Sq  in  such  a  way  that  the  new  Q  wull 
not  be  indefinite,  yet  still  move  the  estimator  eigenvalues  in  the  desired  direction. 

The  constraints  on  the  Q  matrix  being  positive  definite  correspond  to  constraints 
on  the  principal  minor  determinants  of  Q, 

A,  >  0  for  f  =  1,2, ...  ,n  (29) 


w here  A,  is  the  ith  principal  minor  determinant  of  Q.  The  positive  definite  boundaries 
are  the  surfaces  A,  =  0,  which  are  nondinear  constraint  surfaces  in  the  /?n<n+1P2  space. 

Consider  a  change  vector  6q  that  causes  Q  to  cross  a  positive  definite  boundary, 
and  let  the  first  constraint  that  is  violated  be  constraint  k,  i.e.  A*;  =  0  at  this  point. 
The  normal  to  the  tangent  plane  of  the  constraint  surface  at  the  point  of  crossing  is 
found  as  V(Ajt).  'I  his  gradient  vector  is  always  in  a  direction  of  increasing  A*,  that 
is,  back  into  the  positive  definite  region  of  space.  The  elements  of  the  gradient  vector 
are  the  matrix  minors  of  Q  for  the  appropriate  element,  and  can  be  easily  evaluated 
numerically  at  the  point  of  interest  by  successive  determinant  evaluations,  or  for  low 
order  problems  by  analytical  expressions. 

Consider  that  we  are  already  at  a  positive  definite  boundary  At  =  0.  any  further 
moves  to  accomplish  eigen  placement  should  be  done  along  this  boundary.  As  a  linear 
approximation  to  the  boundary  we  move  in  the  tangent  plane,  which  is  completely 
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determined  by  the  gradient  vector  at  the  point  of  contact  with  the  surface.  To  remain 
within  our  minimum  norm  approach,  we  simply  append  another  constraint  to  our 
constraint  equations  (16),  which  then  become 


M  i  Sq 


M 

V(A*)r 


6q1  = 


6  A 
0 


=  6X 


(30) 


and  the  solution  will  be 


6q 1  =  M\6\' 


(31) 


The  new  position  may  still  not  be  in  the  positive  definite  region  due  to  the  con¬ 
vexity  of  the  constraint  surface  at  the  point  of  contact,  so  one  last  move  must  be 
made  back  into  the  positive  definite  region.  The  most  direct  method  would  be  to 
simply  move  in  the  direction  of  the  gradient  vector  at  our  new  point  until  we  cross 
the  constraint  surface.  However,  this  can  lead  to  problems  if  this  move  is  done  with¬ 
out  reference  to  any  eigenplacement  or  minimum  norm  requirements.  In  practice  we 
find  that  in  some  regions  the  constraint  gradient  is  very  nearly  aligned  with  the  direc¬ 
tion  of  greatest  change  of  A.  Thus  even  small  changes  back  to  the  constraint  surface 
sometimes  produce  objectionable  changes  in  system  eigenvalues. 


In  parallel  with  our  algorithm  so  far  therefore,  we  simply  append  a  constraint 
to  (16)  as  before,  but  now  we  wish  no  movement  in  the  eigenvalues  while  increasing 
the  principal  minor  in  question.  Suppose  the  principal  minor  determinant  at  our  new 
point  is  still  some  negative  quantity  —  /?2.  Then  the  problem  becomes 


M\8q2 


(32) 


which  can  be  solved  for  the  second  vector  Sq2.  Finally  then,  the  new  change  vector 
becomes 


6q  =  6qt  +  8q  2 


(33) 
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A  two-dimensional  representation  of  what  is  occurring  at  a  constraint  boundary  is 
given  in  Figure  2. 


Although  it  is  necessary  and  sufficient  that  both  Q  and  P  be  positive  definite  for 
IIe(s )  to  be  strictly  positive  real,  in  practice  one  need  only  consider  the  constraint  on 
Q  because  of  the  way  the  filter  matrix  Kp  is  evaluated.  The  initial  Q  is  chosen  such 
that  P  >  0.  If  P  is  to  become  indefinite,  then  some  small  step  must  cause  P  to  go 
through  the  singular  point  det(P)  =  0.  However,  as  P  comes  closer  to  the  singular 
point,  P_1  will  become  very  large,  hence  I\p  will  become  very  large  due  to  (11).  This 
will  tend  to  significantly  change  the  eigenvalues  of  Pe,  and  the  constraint  that  the 
eigenvalue  placement  be  satisfied  for  that  particular  iteration  tends  to  alter  Q  such 
that  P  decreases,  thus  moving  P  further  from  being  singular. 


3.1.3  Example  1 


We  illustrate  the  procedure  outlined  above  with  a  simple  second  order  example.  Con¬ 
sider  the  system  whose  reduced-order  model  has  been  chosen  as 


Pm  = 


0  1 
2  0  ’ 


Hm  =  [  0  1 


A  feedback  matrix  has  been  chosen  as  K  =  |  0  2 
(Pm  -  OmI<)  ^  —1.0  ±jl.O. 


which  gives  the  eigenvalues  of 


The  estimator  gain  matrix  was  chosen  initially  to  fulfill  the  strict  positive  real 
requirements  of  the  theory  from  a  nominal  choice  of  a  symmetric  Q  as 


where  only  the  upper  triangular  part  of  the  symmetric  Q  is  shown.  This  Q  places  the 
eigenvalues  of  Fe  —  (Pm  —  I\pHm)  at  —0.7628  iy  l.4021.  From  here,  the  eigenvalue 
placement  algorithm  was  applied  to  move  the  eigenvalues  of  Fe  to  — 2.4±j'2.1,  (at  this 
stage  an  arbitrary  eigenplacement).  The  total  8\  desired  is  then  1.6372  ±j 0.6979, 
which,  it  is  easily  checked,  is  well  outside  the  range  where  the  linear  analysis  used  is 
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applicable.  The  step  size  is  an  upper  bound  to  the  allowable  desired  change  vector 
<5A  at  any  step,  and  was  set  to  0.002  for  this  case.  This  step  size  is  small  enough  to 
ensure  that  we  are  within  the  linear  region  of  our  algorithm,  and  can  be  estimated 
by  comparing  the  desired  and  actual  changes  for  one  step. 

After  a  number  of  small  steps,  the  desired  poles  are  reached,  at  which  point 

Q  _  6.S095  5.1317  ' 

V  ~  8.0533 


During  the  steps  to  this  set  of  eigenvalues,  no  positive  definite  constraint  boundary 
was  encountered.  From  this  position,  a  new  set  of  desired  poles  was  specified  as 
—2.75  ±  j2.3,  and  the  eigenvalue  placement  algorithm  applied  again.  Figure  3  shows 
the  resulting  movement  of  the  eigenvalues.  A  positivity  boundary  is  met  during  the 
placement,  but  the  algorithm  compensates  by  appending  the  extra  constraints  to  the 
sensitivity  matrix,  as  in  (30)  and  (32). 


At  the  desired  eigenvalues,  the  final  Q  was 


Q  = 


4.1335  3.7837 
7.4637 


Note  that  for  this  second  order  problem,  once  the  extra  constraint  is  added,  the 
modified  M  becomes  square,  and  at  each  step  the  change  vectors  and  6 A2  become 
uniquely  determined. 

3.1.4  Example  2 

The  second  example  presented  here  is  for  control  of  the  DRAPER  I  tetrahedral  truss 
structure  [11],  shown  in  Figure  4,  which  has  been  used  by  various  authors  as  a  generic 
spacecraft  model.  The  performance  is  measured  as  the  ability  to  control  the  x  and 
y  positions  (line-of-sight  error)  of  the  top  vertex  to  zero.  The  right-angled  bipods 
that  connect  the  truss  to  the  ground  take  on  the  dual  purposes  of  force  actuators  and 
velocity  sensors,  hence  the  collocation  property  is  satisfied  and  the  model  transfer 
function  will  be  positive  real. 
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A  fourth  order  model  of  the  structure  is  formed  using  the  two  modes  of  lowest 
frequency,  which  have  frequencies  =  1.342  rad/s  and  u’2  =  1.665  rad/s  respectively. 
Our  system  matrices  therefore  become,  in  modal  form 

F  _[  [0]  [/] 
m~[-V]  [o] 


where 


M  _ 

^1 

1 - 

0 

r  J  ~ 

0 

w2  J 

There  are  six  actuator/sensor  pairs,  with  Gm  —  H The  Gm  matrix  is  evaluated  from 
the  known  modal  matrix  and  structure  geometry  (see,  for  example,  reference  [5]). 

The  feedback  gain  matrix  K  was  chosen  to  minimize  the  index  of  performance 
given  in  (7)  with  A  =  diag  {103}  and  B  —  I.  The  resulting  poles  of  Fci  =  ( Fm  —  GmK) 
are  at  —1.3177,  —1.4931,  —2.7169  and  —4.5072.  For  the  estimator,  a  set  of  desired 
closed-loop  pole  locations  was  specified  as  —4.0  ±  j2.5  and  —4.5  ±  j3.0.  The  estima¬ 
tor  gain  matrix  is  chosen  initially  using  the  positive  real  equations  from  a  nominal 
diagonal  Q  =  diag  {30.0},  which  places  the  eigenvalues  of  Fe  at  —0.6555  ijl.2321 
and  -1.5480  ±j'0. 8986.  The  eigenvalue  placement  algorithm  was  applied  using  a  step 
size  of  0.001.  The  eigenvalue  paths  are  shown  in  Figure  5,  and  again  we  see  that  we 
meet  a  positive  definite  constraint  during  the  placement. 

At  the  desired  eigenvalue  locations,  the  Q  matrix  is 

23.929  6.1181  29.916  6.3219  ' 

24.373  3.7502  31.057 
38.780  4.0889 
41.758  . 


It  is  important  to  note  that  this  is  not  a  unique  Q  to  achieve  the  desired  eigenplace- 
mont.  Starting  from  a  different  location  in  the  complex  s-plane,  or  even  using  a 
different  step  size,  may  give  a  very  different  final  Q. 
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3.2  General  Non-Linear  Optimization  Method 

For  some  problems,  the  Eigenvalue  Sensitivity  Method  just  described  was  found  to 
exhibit  several  drawbacks.  First,  in  some  cases  the  steps  suggested  to  choose  Sq 
will  not  work.  This  can  occur  for  example,  when  the  desired  eigenvalue  locations 
lie  outside  the  region  where  positive  real  control  is  possible.  Second,  in  some  cases 
investigated,  the  surfaces  of  constant  eigenvalue  and  surfaces  of  constant  principal 
minor  determinant  are  almost  parallel  at  some  points  in  the  iteration.  In  this  case, 
<5A  must  be  chosen  very  small  or  else  the  Sq  calculated  may  be  very  large.  This 
problem  may  sometimes  be  circumvented  by  starting  from  a  different  location  in  the 
s-plane,  but  there  is  no  general  theory  available  to  choose  alternate  starting  points. 
Third,  the  linear  solution  may  also  exhibit  very  slow  convergence  in  regions  where  the 
constraint  surface  has  high  curvature.  These  regions  are  not  known  a  priori  because 
of  the  high  dimensionality  and  complexity  of  the  constraint  surfaces.  In  such  a  case, 
step  sizes  may  need  to  be  very  small  in  order  that  the  linear  approximation  be  valid. 
For  larger  step  sizes,  the  non-linear  terms  in  the  expansion  for  the  constraints  become 
significant. 

To  resolve  some  of  these  difficulties,  an  alternate  form  of  the  pole  placement 
optimization  algorithm  is  proposed.  Rather  than  seek  a  minimum  norm  on  “ Sq ”  for 
a  specified  eigenvalue  change,  we  now  seek  only  to  minimize  the  norm  of  the  error  in 
eigenvalue  placement,  i.e. 


T  j  =  E  I  A,  -  A,  I2  (34) 

1=1 

For  most  problems  this  minimum  is  zero,  and  may  be  achieved  for  a  wide  range  of 
possible  Q  matrices.  To  further  simplify  our  search  for  a  positive  definite  symmetric 
Q,  we  search  rather  for  a  factored  (n  x  n)  matrix  S,  where 

Q  =  STS  (35) 

If  desired,  a  special  form  of  S  (e.g.  a  triangular  5)  instead  of  a  full  5  may  be  used 
to  reduce  the  number  of  variables.  The  definiteness  constraints  are  eliminated  and 
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replaced  by  t ho  sole  constraint 


del  (.S')  ^  I) 


(3(i) 


which  in  fact  will  usually  be  a  non- active  constraint.  It  will  only  become  active  during 
the  solution  proceedure  if  the  set  of  desired  poles  lies  outside  the  positive  real  region, 
or  if  a  traverse  across  the  positive  real  boundary  is  requested. 

As  in  the  previous  method,  only  those  eigenvalues  required  to  be  placed  need  be 
included  in  the  sum  in  (34),  and  if  the  imaginary  parts  of  the  eigenvalues  are  free, 
then  these  components  can  be  deleted  from  the  objective.  This  objective  has  the 
advantage  over  that  used  previously  in  that  if  the  set  of  desired  eigenvalues  is  not 
a  subset  of  the  achievable  set,  then  the  minimization  remains  well  posed  and  will 
converge  to  a  set  of  poles  as  close  as  possible  to  the  desired  set.  In  principle,  this 
technique  could  be  used  to  numerically  evaluate  the  positive  real  region  for  particular 
pole  groupings. 


Derivatives  of  the  performance  index  (34)  are  no  more  difficult  to  evaluate  than 
the  eigenvalue  sensitivities  evaluated  previously.  Using  the  same  methodology,  the 
first  derivative  of  the  objective  function  (34)  with  respect  to  the  elements  of  5  is 


OJ 

dSjk 


2£|  A.-A, 


f  d\j 
OS 


jk, 


(37) 


where,  in  parallel  with  the  first  algorithm. 


o>A.  __  (  dP  \ 

—  V ,  P  I  ^  ^  I  A  fr 


dSjk 


dS 


jk  , 


P  is  the  positive  definite  solution  to 


PFci  +FjlP  =  ~  [St5  -  I<rHm  -  HlK 


(38) 


(39) 
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dP 

osjk 

is  a  matrix  obtained  by  differentiating  equation  (39),  which  gives  the  Lyapunov  equa¬ 
tion 


dP 

dS]k 


Fa  +  Fl 


dS 

dSjk 


T 

S  +  ST 


(40) 


The  solution  to  the  above  constrained  optimization  problem  was  found  using  the 
subroutine  package  GRG2  (Generalized  Reduced  Gradient)  developed  by  Lasdon  and 
Waren  [12]  in  1982. 

3.2.1  Examples  1  and  2  Revisited 

When  examples  1  and  2  are  solved  again  using  this  second  formulation,  the  poles  are 
again  placed  as  desired,  and  the  final  Q  matrices  obtained  are  as  follows. 


5.5826  4.8579 
8.1710 

'  24.131  -3.9992  31.307  -0.7582  ' 

83.110  -8.7024  90.621 

41.045  -5.1034 
100.69 


Note  that  these  Q  matrices  are  very  different  from  those  evaluated  using  the 
first  algorithm,  yet  accomplish  the  same  pole  placement.  This  underscores  the  non¬ 
uniqueness  of  the  solution  matrix  Q,  as  well  as  the  non-uniqueness  of  Kp  (in  general). 
Since  for  this  case  we  used  the  GRG2  code,  no  external  search  path,  such  as  that  in 
Figures  4  and  5,  was  found. 


For  the  two  methods  presented,  neither  appears  to  produce  a  consistent  control 
system  performance  superiority,  although  numerically  the  latter  approach  is  clearly 
more  efficient  and  robust. 
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3.2.2  Example  3 


The  last  example  is  a  sixth  order  model  ol  t.he  DRAPER  I  truss  structure.  This 
example  was  used  to  generate  the  eigenvalue  locus  given  in  Figure  1  which  illustrates 
the  estimator  eigenvalue  placement  when  using  a  diagonal  Q  —  ql .  In  this  example 
we  illustrate  the  aditional  flexibility  in  eigenvalue  placement  afforded  by  allowing  an 
arbitrary  (non-diagonal)  Q  =  S1  S.  The  modes  kept  in  uhe  model  were  modes  1, 
2  and  5  (in  [3]  these  modes  were  determined  as  “most”  important).  Mode  5  has 
frequency  u-'s  =  3.398  rad/s.  The  feedback  gain  matrix  I\  is  chosen  to  minimize 
(7)  with  A  —  diag  { 1 03 }  and  B  =  /.  The  resulting  poles  of  Fci  are  at  —1.2754, 
-1.4931,  -2.7169,  —5.8560.  and  —0.6509  ±  j 3.0028.  The  estimator  gain  matrix  is 
chosen  from  an  initial  choice  of  5  =  So  =  diag{>/50},  which  places  the  eigenvalues  of 
/•;  at  -0.347  ±  j  1.322,  -0.296  ±  j  1.781.  and  -1.289  ±  j3.396. 

A  set  of  desired  pole  locations  is  specified  as  —1.70  ±  j  1 .70,  — 1.80  ±  jl.SO.  and 
—  3.00  ±  j 3.50.  This  set  of  desired  poles  is  chosen  rather  arbitrarily,  the  point  of 
the  exercise  being  to  show  the  efficiency  of  the  eigenvalue  placement  algorithm.  In 
previous  work  [3],  we  had  found  that  the  inability  to  do  this  with  a  diagonal  Q  (as 
shown  in  Figure  1)  led  to  poor  controller  performance. 

(  sing  the  second  formulation,  this  problem  was  solved  using  the  optimization 
package  GRG2,  and  at  the  desired  eigenvalue  locations,  the  Q  matrix  is 


-19.90 

-4.272 

29.70 

-21.36  ‘ 

-11.86 

27.70 

33.84 

-0.953 

109.3 

3.445 

-66.26 

52.38 

22.57 

9.321 

-2.100 

60.10 

-35.29 

30.38 

Figure  6  shows  the  eigenplacement  with  respect  to  the  initial  eigenvalues  of  the  estima¬ 
tor.  It  is  obvious  that  a  significant  increase  in  the  damping  ratios  for  the  eigenvalues 
of  F,  has  been  achieved  while  leaving  the  damped  frequencies  almost  constant.  This 
solution  required  approximately  60  CPU  seconds  on  an  Amdahl  470  computer. 
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The  question  remains  whether  or  not  this  eigenvalue  placement  results  in  improved 
system  performance  over  a  diagonal  Q.  To  address  this  question,  we  form  a  positiv? 
real  controller  for  this  system,  using  a  diagonal  Q  of  the  form  Q  =  ql ,  which  has 
the  same  “control  cost”  as  the  controller  using  the  general  pole  placement  presented 
here.  The  control  cost  here  is  defined  as 

CC  -  ( utBu )  dt  =  xlQZxaa  (41) 

where  xj  =  j  xT  xV  j,  Z  satisfies 

ZFa  +  FjZ  =  —Ka  (42) 


and 


F  -GI< 
a  ~  I<fH  Ft 


(43) 


‘  0  0 
0  I\tRK 


(44) 


Tlie  cost  is  evaluated  for  a  particular  set  of  initial  conditions  xao.  In  this  case.  xao 
was  chosen  as  a  unit  displacement  in  the  x-y  plane  of  the  top  vertex  of  the  truss, 
equally  shared  in  the  x  and  y  directions. 

If  there  is  no  model  error,  i.e.  if  F,  G  and  H  are  Fm,  Gm  and  Hm  respectively, 
then  the  predicted  eigenvalues  of  Fci  and  Fe  are  indeed  the  true  system  closed-loop 
eigenvalues.  For  this  case,  shown  in  Figure  7,  the  system  whose  poles  have  been 
placed  to  improve  their  damping  exhibits  much  better  error  response  that  the  system 
whose  poles  were  chosen  on  the  basis  of  a  diagonal  Q,  even  though  the  same  control 
energy  is  used. 

If,  however,  we  use  the  full  system  matrices  for  F ,  G  and  H,  we  introduce  ad¬ 
ditional  unmodeled  modes,  as  well  as  control  and  observation  spillover.  For  this 
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case,  we  observe  the  time  response  as  shown  in  Figure  8.  The  significant  increase 
in  the  damping  ratios  of  the  estimator  eigenvalues  does  not  translate,  in  this  case 
at  least,  to  a  significant  improvement  in  the  LOS  error  response.  For  this  problem, 
the  residual  modes  in  the  system  (those  modes  not  in  the  controller)  are  excited  by 
this  disturbance,  anti  the  improved  damping  in  the  estimator  modes  does  not  affect 
them.  This  is  a  manifestation  of  the  general  problem  of  modal  approximation  and 
model  order  reduction.  When  estimator  feedback  is  used  as  the  basis  for  controller 
design,  the  predicted  closed-loop  poles  are  simply  the  poles  of  Fci  and  Fc.  In  practice 
when  a  reduced-order  model  is  used  to  design  the  estimator  as  in  our  examples,  the 
actual  closed- loop  poles  may  be  quite  unpredictable,  depending  on  the  fidelity  of  the 
approximate  model. 

4  Conclusions 

A  new  technique  has  been  presented  for  the  determination  of  positive  real  transfer 
matrices  with  desired  eigenvalues.  The  technique  casts  the  problem  as  a  constrained 
optimization  problem  which  is  solved  iteratively  using  gradient  methods.  Two  solu¬ 
tion  algorithms  have  been  presented  which  differ  in  the  optimization  problem  formu¬ 
lation.  The  second  algorithm,  which  is  based  on  nonlinear  optimization,  was  found  to 
be  superior  for  the  examples  used  in  terms  of  computer  time,  overall  convergence  char¬ 
acteristics  and  robustness.  The  new  algorithms  were  successfully  applied  to  a  simple 
second  order  problem,  and  to  fourth  and  sixth  order  models  of  the  DRAPER  I  struc¬ 
ture.  Significant  changes  in  the  estimator  eigenvalue  positions  have  been  achieved 
while  retaining  the  positive  realness  of  the  transfer  matrix.  This  technique  promises 
to  give  greater  flexibility  in  design  when  using  positive  real  control. 
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Figure  3:  Convergence  History  for  Estimator  Eigenvalues  for  2nd  order  model 
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Figure  6:  Initial  and  Final  Eigenvalue  Placement  for  6th  order  model  of  the  DRAPER  I 
Structure 
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Figure  7:  LOS  Error  Response  Using  a  Reduced-Order  System  Approximation 
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F'guie  b.  LOS  Error  Response  ot  Full  System 
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